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The quantum XX chain - or rather ring - is studied as a toy model of an interface. Two transverse 
field patterns are used to define the interface, on the one hand a staggered field, on the other hand 
a step-like configuration, from —h to +h. The interface leads to Friedel oscillations and proximity 
effects, in particular close to the quantum phase transition of the bulk, which is a metal-insulator 
transition in the fermionic language. The most prominent interface effects appear for odd-numbered 
rings, for which - in contrast to even chains - the ground state is doubly degenerate. In the regime 
where the bulk energy spectrum is gapped a level appears close to midgap, with a wave function 
localized in the region of the interface. The two members of the ground state doublet have two 
different particle number parities and spin components S z = ±|. They also have different energy 
levels and thus the degeneracy does not originate from different occupancies of a rigid band structure, 
but rather from a global symmetry. The unitary transformation linking the two degenerate ground 
states resembles a Majorana operator. Coherent superpositions of the two states may be suitable 
candidates for well protected qubits. 

PACS numbers: 75.10.Pq 



I. INTRODUCTION 



Spin chains are favorite models for studying fun- 
damental quantum phenomena. Thus the XY chain, 
where only the x and y components of the spin opera- 
tors are coupled, has been used for discussing quantum 
phase transitions ^^, quantum dynamics^^ and quan- 
tum entanglemenlPA It has also served as an illustra- 
tive case displaying the effects of randomness on critical 
properties^ and as a handy example for scrutinizing vari- 
ational ground states^. 

In this paper we use the model for discussing charac- 
teristic effects produced by an interface. A multitude of 
interesting phenomena can occur at the interface sepa- 
rating two different materials, such as proximity effects 
- the spillover of the order present in one material into 
the other within a certain spatial region - Friedel oscil- 
lations - a periodic variation of electronic charge density 
in a metal close to the interface with an insulator - or 
the generation of a two-dimensional electron gas at the 
interface between a doped and an undoped semiconduc- 
tor, a prerequisite for the Quantum Hall Effect. Recently 
metallic - and even superconducting - states generated at 
the interface between a band and a Mott insulator have 
stimulated a whole new field of research^. Similarly, 
in topological insulators conducting states are produced 
at an interface (or at the surface), with rather exotic 
properties^. 

We model the interface by choosing particular config- 
urations of the transverse fields, as illustrated in Fig. [T] 
For the staggered field an interface is created by joining 
the ends of an odd-numbered chain. In the step model 
the interface is produced by applying opposite fields in 
the two halves of the chain. 

We consider the isotropic version of the model, the XX 
chain - or rather a closed ring - coupled to a transverse 
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FIG. 1: (Color online) Examples of inhomogeneous field configu- 
rations: a) staggered field, b) step model with a sharp interface. 



field. The Hamiltonian is 

L 

where S", a — x,y,z, are spin | operators, S" +1 = Sf 
and J is chosen positive. The Jordan- Wigner transfor- 
mation 
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to spinless fermion operators Cj , Cj leads to a Hamiltonian 
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of non-interacting particles 
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where the particle number N — j c ] c 3 i s related to the 



z component of the total spin as 
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For odd N, the Hamiltonian ^ is a simple tight-binding 
model with an inhomogeneous chemical potential. For 
even N, the hopping term between sites L and 1 has a 
different sign than the others. For very large chains this 
should have negligible effects on bulk properties. How- 
ever, as we will see later, this boundary term is essential 
for the ground state degeneracy in odd-numbered chains. 

Before discussing interface effects we recall the essen- 
tial ground-state properties in the absence of an inter- 
face. Specifically, we consider the two cases of constant 
and alternating fields (for an even number of sites in the 
latter case). Some mathematical details are presented in 
Appendix [X] 



For a constant field, hi 



, L, the main 



properties of the model have been widely explored, ever 
since the seminal work of Lieb, Schultz and Mattis^. 
The single-particle spectrum J cos k — ft is negative for 
all values of the wave vectors, — 7r < k < 7r, if h exceeds 
J; for h < — J, the spectrum is positive. For h = the 
ground state corresponds to a half-filled band with Fermi 
wave vectors ±7r/2. The filling changes as a function of 
ft, with Fermi wave vectors moving inwards for ft > and 
outwards for ft, < 0: they disappear at ft, = ± J, where a 
Lifshitz transition^ from a partially filled band to either 
a completely filled or a completely empty band occurs. 
The magnetization per site 
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< ft < J, 
h> J 
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(and correspondingly for ft, < 0) is shown in Fig. [2j to- 
gether with the susceptibility x = dm/dh. The two crit- 
ical points h — ±J mark the transitions from a partial 
magnetization to fully aligned moments. 

For a staggered field, hj = (— f) 3 ft,j = 1, —,L, the 
single-particle spectrum has a gap 2\h\ separating a va- 
lence from a conduction band for an arbitrarily small field 
amplitude ft. Moreover there is long-range antiferromag- 
netic order with local moments given by 
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where K,(x) is a complete elliptic integral of the first kind 
and A = ft/ J. Fig. [3] shows both the order parameter 




FIG. 2: (Color online) Magnetization (lower line) and susceptibil- 
ity (upper line) for a homogeneous transverse field. 



(the staggered magnetization (— l) J mj) and the suscep- 
tibility x = ( — l) 3 dnij/dh. The order parameter grows 
first linearly as a function of the field and saturates for 
| ft | 3> J. The susceptibility diverges at the critical point 
ft = 0. 

The ground state can also be readily obtained if both 
homogeneous and staggered fields are included. Besides 
phases with fully aligned or alternating spins a third 
phase is found where both orders coexisfr^l 




FIG. 3: (Color online) Staggered magnetization (— l) J mj, pass- 
ing through at h = and saturating at ±^ for \h\ — > oo, and 
susceptibility Xi diverging for h — > 0. 



From now on we will mainly concentrate on interfaces, 
such as the step model or odd-numbered rings with stag- 
gered transverse fields. Sometimes more general inhomo- 
geneous field patterns will be considered, subject to the 
constraint 



(7) 
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This condition of vanishing average field will also be ap- 
plied to the two particular configurations of Fig. [TJ where 
it holds automatically for even-numbered chains. To sat- 
isfy the constraint for an odd number of sites we assume 
the field to vanish at a single site (at the site of the in- 
terface for the step model), 
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Instead of inhomogcneous fields one could also consider 
inhomogeneous exchange constants Jj. Their effects are 
however expected to be much weaker than those of an 
inhomogeneous field. Consider an even-numbered chain 
with exchange constants varying arbitrarily in sign, but 
not in absolute value. Applying the canonical transfor- 
mation S? -> —Sf, Sj -> -Sj and Sj -> 5| at appro- 
priate places then leads to a model with homogeneous 
exchange. In contrast, one cannot simply gauge away 
varying signs of transverse fields without modifying other 
terms in the Hamiltonian. 

The paper is organized as follows. Symmetry consider- 
ations are presented in Section II, which are very useful 
for the XX chain with interface, despite the fact that 
translational invariance is explicitly broken. For general 
fields vanishing on average the concept of particle num- 
ber parity can be used to determine the ground state 
degeneracies. For an even number of sites and more re- 
stricted field patterns, namely those that are odd under 
reflection about the chain center, we find a particular 
SU(2) symmetry involving operators O x ,O y ,O z . While 
O z is essentially equal to the z component of the total 
spin, the other two operators generate mappings con- 
sisting of both a reflection and an electron-hole trans- 
formation. Section III discusses the step model for an 
even-numbered chain. The critical point separating a 
"metallic" from an "insulating" phase is at \h\ = J, in 
close analogy to the Lifshitz transition for a homogeneous 
field. For \h\ < J, pronounced Friedel oscillations in the 
local moments (or the particle density) occur. Both their 
characteristic wavelength and their intensity depend on 
the field strength. Chapter IV deals with the intriguing 
case of odd-numbered chains, both for staggered fields 
and for the step model. In both cases an isolated level 
appears in the gap, which is not exactly at zero energy. 
Nevertheless the ground state energy does not depend on 
the occupancy of this level, in perfect agreement with the 
symmetry arguments of Section II. Although the SU(2) 
symmetry of even chains no longer holds, the unitary 
transformations introduced in Section II can be general- 
ized to odd chains. One of these transformations links 
the two degenerate ground states and has the form of a 
Majorana operator. A coherent superposition of these 
two states, a possible qubit, has a magnetic moment in 
the plane, with some angle and an amplitude < i. Sec- 
tion V presents a brief summary as well as a discussion 
of possible extensions, such as the XXZ model or the XY 
chain with anisotropic exchange. 



II. SYMMETRY CONSIDERATIONS 

A. Particle-number parity 

The Hamiltonian (13J) depends on the particle-number 
parity P — exp(i7r7V)7a concept that has been success- 
fully applied in other contexts, such as even-odd effects 
in small metallic grains^. This operator commutes with 
the Hamiltonian and therefore both the single-particle 



levels and the many-particle eigenstates can be grouped 
in two parity sectors according to the eigenvalues of P 
(±1). For transverse fields satisfying the constraint ^ 
the Hamiltonian has also particle-hole symmetry in the 
following sense. The transformation cj — > cj, cj — > Cj, 
j = 1, L, simply changes the sign of the first two terms 
of the Hamiltonian. As to the boundary term, one has to 
take into account that at the same time N is transformed 
to L — N and therefore P — s- P for even L and P — > — P 
for odd L. Therefore this particle-hole transformation 
amounts to the mapping 

(H P\ ^ i (~ H > P }> L even ' (R\ 
(ff,P)->| (- H ,-P), L odd. (8) 

Because this is a canonical transformation, the spectrum 
does not change. For even chains the energy levels oc- 
cur in pairs {e v , — e u ) in both parity sectors. For odd 
chains the energy level e v of one parity sector has its 
partner —e v in the opposite parity sector. To determine 
the ground state one has to explore both sectors. Setting 
aside the ambiguity of zero-energy levels, one quickly re- 
alizes that the rule of simply occupying all the negative- 
energy levels is not always consistent. In fact, the calcu- 
lation of the single-particle spectrum for a given parity 
sector may produce a number of negative-energy levels 
that disagrees with the parity chosen initially. For even 
chains the number of negative-energy levels is equal to 
L/2, in both parity sectors, and one expects the ground 
state to have parity P = exp (in L/2), for which it is con- 
sistent to occupy the negative-energy levels and to leave 
empty the positive-energy levels. Thus the ground state 
is expected to be unique for even chains. 

For odd chains the situation is slightly more compli- 
cated. If there are N negative-energy levels in one sector, 
there are L — N negative energy levels in the other, but 
this number N may or may not agree with parity, as will 
become more clear for the explicit examples of Section 
IV. Let us first discuss the case where parity and num- 
ber of negative-energy levels are consistent. The ground 
state energy in the + parity sector is then given by 

£«= £ 4+). (9) 

vd +) <0 

Together with the relation J2„ £ ^ = 0, a simple conse- 
quence of the constraint ^ ■ hj — 0, and the fact that to 

each energy level there exists a level ei ' = — e! + \ 
we obtain 

e {+) = - E 4 +) = E ^ = &-'>■ (io) 

Therefore there are two orthogonal states with the same 
energy. If the particle-number parity and the number of 
negative-energy levels are not consistent, one has either 
to add a particle to the lowest unoccupied positive-energy 
level or to remove a particle from the highest occupied 
negative-energy level. Particle-hole symmetry together 
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with the vanishing trace of eigenvalues can then again be 
used to show that the many-particle states in the two sec- 
tors have the same energy. We conclude that the ground 
state for an odd-numbered chain is degenerate. 



B. SU(2) symmetry (L even) 

Consider an even-numbered chain, L = 2M, with fields 
that are odd under reflection about the chain center, i.e. 



hi = —h 



L+l-j, 



j = l,...,M. 



(11) 



The two field patterns of Fig. [T] satisfy this relation. We 
now introduce the operators 
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They satisfy the commutation relations of an angular mo- 
mentum 

[0 X , Oy] = iO Z , [Oy, Z ] = iO X , [0 Z ,0 X ] = iOy (13) 

and commute with the Hamiltonian ^ for field config- 
urations satisfying Eq. (11). Therefore our model has 
SU(2) symmetry. We note that the XY Hamiltonian ([3| 
has of course not the SU(2) symmetry with respect to the 
true spin operators. Rather the operators ( 12 ) are similar 



to the 77 spin introduced by Yang and Zhang in the case of 
the Hubbard modeP^. The possible eigenvalues of O z are 
integers between — M and +M. Thus the many-particle 
states may have rather large degeneracies, in contrast to 
the single-particle states, which are expected to be non- 
degenerate. Certain (many-particle) eigenstates can then 
be constructed in a similar way as those of the ordinary 
angular momentum operator. For instance, starting with 
the eigenstate for N = L, where all sites are occupied, 
we can generate an eigenstate for N = L — 2 by applying 
the lowering operator 
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0-=O x 
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In a similar way one finds an eigenstate for N = 2 by 
applying 

M 

0+ = O x + iO y = V2J2(-^4 c I+i-j ( 15 ) 

3=1 

to the ground state for N — 0, where all sites are unoc- 
cupied. 



C. Symmetry transformation 

The operators introduced above can be used as gen- 
erators of symmetry transformations. As a particular 
example we consider the unitary operator 



U : 



e si 



(16) 



which induces a particle-hole transformation plus reflec- 
tion, 



Uc 3 U^ = it-lM+i-j 

ucL+t-jrf = i{-iy c ), 



(17) 



where j = 1, ...,M and L — 2M. It leaves the Hamilto- 
nian ([3| invariant if Eq. (Ill holds. One also finds the 
simple relation 



M 



tf 2 = IJ[2(n J --n L+1 _ i ) a -l] 
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where 



- J, 



This operator is both hcrmitian and 



unitary and therefore its possible eigenvalues are ±1. 
As an example we choose a state with an even num- 
ber of particles distributed over the L sites such that 
n-j = UL+i-j, j = 1,...,M. In this case U 2 = (-1) M . It 
is easy to check that any permutation of particles does 
not change this value and therefore any state of an even 
number of particles (or of even particle number parity) is 
an eigenstate of U 2 with eigenvalue (— 1) M . Similarly, for 
an odd number of particles we can start with a configura- 
tion where the relation rij = riL+i-j, j — X,...,M holds 
for all sites except one. Therefore one of the factors in 
Eq. (18) is equal to 1 whereas all the others are -1. The 
eigenvalue of U 2 for any state with odd number parity 
is — ( — 1) M . Thus any iV-particle state of a chain with 
L = 2M is an eigenstate of U 2 with eigenvalue (—1) N+M . 
We conclude that up to a phase factor of (— 1) M the op- 
erator U 2 is identical to the particle number parity P. 



III. FRIEDEL OSCILLATIONS (L = 2Af) 

The step model of Fig. [T]can be viewed as the union of 
two chains with different homogeneous transverse fields. 
If interface effects are neglected, each part will undergo 
separately a quantum phase transition, one at the left 
critical point of Fig. [2] the other at the right point. 
This can be quantified as follows. We introduce fermionic 
operators Qj , bj through 



aj , j = l...,M 

h, j = M + l, 



(19) 



and write 
H 
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where H' consists of boundary terms proportional to J 
and om+i = a\, &a/+i — b±. If we neglect H' we find two 
independent chains with homogeneous fields —h and +h, 
respectively, and with single-particle spectra Jcosfc ± h. 
For \h\ > J one band is empty the other one is full. 
In this limit, the magnetic moments m,j are —\ for j — 
1, M and +| for j = M + 1, L. 



e/J 




in the figure. A closer look shows that for some spe- 
cific field strengths the oscillations are practically ab- 
sent, while they are particularly pronounced for others. 
We attribute this effect to interference between Friedel 
oscillations originating from the two interfaces located, 
respectively, in the middle and at the end of the chain. 

m j 
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FIG. 4: (Color online) Energy levels of the step model for L = 18. 
An energy gap opens for h/.J > 1, signaling the metal-insulator 
transition. 

The boundary terms in H' are negligible for \h\ 3> J, 
but not for \h\ < J. To study the "metallic" region 
we have numerically diagonalized the Hamiltonian. The 
single-particle spectrum is illustrated in Fig. [4] A gap 
exists for \h\ > J, while the system is gapless for \h\ < J, 
as expected. The local moment at site j is given by 

1 
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f,£<0 



*J" I 
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where the matrix elements of the unitary trans- 

formation that diagonalizes the Hamiltonian (y numbers 
the energy levels). The numerical results are illustrated 
in Fig. [5] for several values of the field strength. Only 
the right half of the chain is represented, the left seg- 
ment being simply obtained by inversion with respect to 
the origin of the axes. Pronounced oscillations are ob- 
served with amplitudes that first increase as h increases 
and then quickly disappear at the approach of the criti- 
cal field strength, where the magnetization saturates at 



in j 



+ | in the right half of the chain (and 



■5 to the 



left). Both for h <C J and for h ^> J the magnetic mo- 
ments can be calculated using perturbation theory, as 
shown in Appendix [B] In the appropriate limits the re- 
sults agree well with those shown in Fig. [5] Thus for 
/i > J the expansion in powers of J/h shows that the 
mutual influence of the two different orderings is only 
effective in close proximity to the interface. 

To make contact with the conventional picture of 
Friedel oscillations, we split the Hamiltonian into two 
parts as in Eq. ( 20 1 and consider field strengths smaller 
than J. The first term represents two metallic systems 
with Fermi wave vectors ± arccos (hj J). The surface 
term H' is then expected to produce Friedel oscillations 
with period n/ arccos (hj J). This prediction is in per- 
fect agreement with the periods of the oscillations seen 



FIG. 5: (Color online) Magnetic moments mj of the interface 
model for L = 102 and j = 51, 102. The different lines corre- 
spond to different field strengths; these are, from bottom to top, 
h = 0.00, 0.10, 0.90, 0.95, 0.97, 0.99, 0.993, 0.999. 



IV. ODD-NUMBERED CHAINS 



We turn now to the discussion of odd-numbered chains. 
To preserve the condition Y]j hj — 0, we start from an 
even-numbered chain and add a site at which the mag- 
netic field vanishes. The choice of the site is arbitrary for 
the staggered field while for the step model we insert it 
at the interface. Choosing in both cases the additional 
site at £ = (L + l)/2, we can write the field term of the 
Hamiltonian as 



l-X 



(22) 



where Oj = (— l) 3 for a staggered field and Oj = 1 for the 
step model. 



A. Midgap level 

To understand qualitatively the spectrum, we consider 
first the limit \h\ ^> J and bypass the new site by re- 
placing the hopping terms involving the site I by a single 
hopping between I — 1 and I + 1. We obtain two inde- 
pendent systems, an even chain - with either a staggered 
field or a field step - plus a single site without field. In 
this limit the spectrum consists of a filled valence band, 
an empty conduction band and a single midgap level at 
s = 0. 

The numerical diagonalization of the Hamiltonian con- 
firms this picture. Fig. [6]shows the result for a chain with 
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11 sites and a staggered field at all but one site. The 
spectrum clearly exhibits the two bands plus an isolated 
level, which tends very quickly to midgap as h increases. 
The two figures correspond to two different particle num- 
bers. N is odd in the upper case and even in the lower 
case. This is at odds with the number of negative-energy 
levels, which is even in the upper case and odd in the 
lower case. Therefore in the ground state the midgap 
level is not occupied in the upper case, although its en- 
ergy is slightly negative, while it is occupied in the lower 
case, although its energy is positive. The particle-hole 
symmetry between the two figures is obvious. It implies, 
as explained in Section II, that the total energies of the 
two distinct states are equal, hence the ground state is 
doubly degenerate. For the step model similar results are 
obtained for \h\ > J. 

In the presence of an energy gap of size 2A the wave 
function UjQ of the midgap state is localized, i.e. it 
decreases as exp (\j — with a localization length 

£ fa J I A. For the step model (in the limit L — > oo) 
£ diverges at the critical point. 



b) 




FIG. 6: (Color online) Single-particle spectrum of the odd- 
numbered chain L = 11 with a staggered field. Both parity sectors 
are shown, in a) N is odd while in b) N is even. Note that all 
energy levels are non-degenerate; for one of the outermost levels 
the splitting is barely visible. 



Hamiltonian, but the modified unitary transformation 

U' = e- i ^°*{c i + 4) (23) 

together with a local gauge transformation does. To see 
this, we first notice that the relations (17 1 remain valid 
for j = 1, . . . , I — 1. Moreover we find 



U'dU 



't _ J 



(24) 



The application of this transformation together with the 
replacement eg — > i(— 1) q indeed leaves the Hamiltonian 
invariant. 



In contrast to the operator U of Eq. (16 1, which is 



bosonic, the operator U' is fcrmionic. Moreover it satis- 
fies the relations 



17* = U , (-l) M+N , U' 2 = (-1) 



M+N 



(25) 



and thus is similar to the creation (or annihilation) op- 
erator of a Majorana fermion. 



C. Degenerate ground state 

We show now that the members of the ground state 
doublet are transformed into each other by the unitary 
operator U'. We use the notation |±) for the two states 
which are eigenstates of S z with eigenvalues ±1/2. In 
fact, in view of Eq. Q we get 



S'\±) = (N-^)\±) = ±±\±) 



(26) 



Because U' commutes with the Hamiltonian, the state 
U ,J< \+) has the same energy eigenvalue as the state |+) 
and must be a superposition of |+) and |— ). Moreover 
the relation U'S Z U^ = -S z implies 



S z U r 



-U'^S Z \ 



(27) 



and therefore we can identify with |— ). We con- 

clude that the two degenerate states |±) are transformed 
into each other as 



U' 



u'\ 



(28) 



An arbitrary superposition of these two states is also 
an eigenstate of the Hamiltonian. We may define a qubit 
as the state 



1 



l+>+e l/3 |-» 



(29) 



B. Symmetry transformation 



The operators ( 12 1 can also be used for odd- numbered 
chains, with M — £—1. They no longer commute with the 



The expectation value of S z for this state vanishes, but 
not that of S x or S v . To estimate (P\S a \P), a = x,y, 
we use the limit \h\ 3> J, for which we can cut the con- 
nections between the added site and the rest of the chain 
(the neglected terms are of order J) . We obtain the fac- 
torization 



|$>®|1>, 



|$>® |0) 



(30) 
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where |1) and |0) are the two possible states of the added 
site with occupation numbers 1 and 0, respectively, and 
|$) is the ground state for the Hamiltonian involving all 
the other sites. In this approximation |$) does not de- 
pend on the occupation of the added site. It is convenient 
to split the spin operators as follows 



ga = + go 

where I is the added site and 



S' a =j2 s f- 



(31) 



(32) 



Because of Eq. (26) |$) can be assumed to be an eigen- 
state of S' z with eigenvalue zero. It follows that also the 
expectation values of S' x and S' v with respect to |3>) van- 
ish. It is then easy to see that for the qubit state \/3) we 
obtain the simple result 

(0\S x \p} = ^ COS p, </3|S»|/3) = ± sin/3. (33) 

In the limit \h\ 3> J the double degeneracy originates 
from the existence of a zero-energy level, which may or 
may not be occupied. For smaller fields the degeneracy is 
not simply related to different possible occupancies of the 
zero-energy level, but relies on global symmetries. This is 
quite different from the famous midgap states associated 
with kink solitons in the SSH model 17 , where the zero- 
energy level can be empty, singly or doubly occupied, and 
where neither the midgap state nor the band states are 
affected by a change in occupancy. In our case a change 
in the particle number by 1 modifies the Hamiltonian, 
and therefore both the midgap state and the band states 
are altered. 

It is remarkable that the degeneracy of the ground 
state does not depend on the field configuration except 
that the trace is required to vanish, Y]- hj — 0. If this 
condition is fulfilled, field fluctuations will not lead to de- 
coherence. Moreover, for the two examples studied here 
the doubly degenerate ground state is separated from the 
excited states by a finite gap, provided that the field 
strength is large enough in the case of the step model 
(| h | > J). Thus the conditions for a useful qubit, low 
decoherence and excitation gap, seem to be well satisfied 
by our model. Both the relation ( 26 ) and the subsequent 



discussion remain valid, except that now there is no guar- 
antee for the existence of an energy gap. Therefore the 
state ( 29 ) cannot be used indiscriminately as a well pro- 



tected qubit. 



D. Generalization 

We can easily extend these considerations to more 
general field patterns, namely to those defined by Eq. 
(Ill for j = 1, ..,,£ — 1 together with the condition 
hi = 0. For simplicity we use open boundary condi- 
tions, for which the parity-dependent term in the Hamil- 
tonian Q is absent. In this case the transformation 



Cj -H> ( — lpCL+i-j, j = 1, ...,£, maps H into —H, and 
the single-particle spectrum is symmetric, i.e. energy 
levels occur in pairs (e a , —s a ), except for one level which 
is mapped onto itself and therefore has zero energy. In 
one of the ground states the zero-energy level is empty, 
in the other it is occupied; they correspond to the states 
|+) and |— ) described above. 



V. DISCUSSION 

In this paper we have studied the ground state of the 
XX chain with an interface defined by specific patterns 
of the transverse field. Two types have been investigated 
in detail, a staggered field where an interface is produced 
automatically for odd-numbered rings and a step model 
where the field jumps suddenly from a negative value 
-ft to a positive value +h. The step model illustrates 
nicely the effects of an interface separating two halves 
with different magnetizations, or with different particle 
densities in the fermionic language. In the "metallic" 
regime, \h\ < J, we have found Friedel oscillations which 
become quite pronounced with increasing field strength 
but then quickly die out as the critical point \h\ = J is 
approached. For larger fields a proximity effect remains, 
limited to the sites adjacent to the interface. To realize 
the step model, or at least a smooth step satisfying the 
condition ( 11 ), one could imagine two identical U-shaped 



magnets placed symmetrically with respect to the center, 
but with north and south poles interchanged. 

Odd-numbered rings exhibit some rather intriguing 
phenomena. For any field pattern satisfying the rela- 
tion J2j hj — the ground state is doubly degenerate. 
In the insulating phase of our two interface models this 
degeneracy is linked to a level close to the center of the 
energy gap. For finite lengths the precise location of this 
level shifts as the particle-number parity P = exp (iirN) 
is changed, due to a parity-dependent boundary term. 
Thus the degeneracy for an odd number of sites origi- 
nates from a global symmetry and not simply from the 
occupancy of a zero-energy level. An explicit unitary 
transformation U' has been constructed, which maps one 
ground state to the other. The operator U' has the prop- 
erties of a creation (or annihilation) operator of a Majo- 
rana fermion, and U' 2 is essentially equal to the parity P. 
The doubly degenerate ground state has spin S z = ±| 
and may be a promising candidate for a well protected 
qubit. 

An interesting question is to what extent the present 
results survive if the coupling between z components of 
the spins is also included. In the fermionic representa- 
tion it amounts to considering an interaction term be- 
tween particles on nearest-neighbor sites. This problem 
has been studied numerically both for homogeneous and 
staggered fields^, but to our knowledge not in the pres- 
ence of an interface. The coupling between z components 
breaks the SU(2) symmetry established in Section II, yet 
it does not spoil the invariance of the Hamiltonian un- 



der the transformations ( 16 1 and ( 23 ) for even and odd 
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chains, respectively. Therefore the eigenstates of odd- 
numbered chains remain degenerate. 

A further possible generalization would be anisotropic 
exchange, J x ^ J y . In this case the Jordan- Wigner trans- 
formation leads to a BCS-type model and, as shown by 
Kitaev-^, to the possibility of Majorana fermions at the 
two ends of a chain in a (homogeneous) transverse field. 
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Appendix A: Solution of the XX chain for 
homogeneous and staggered transverse fields 

1. Homogeneous field 

For odd N and a constant transverse field, hj = h,j = 
1,...,L, the Hamiltonian ^ is diagonalized by Fourier 
transformation 



-y 



e ikj c k , 



where k = ^v, -L/2 <v< L/2. We find 



H = ^(J cos k-h)c\ 



Ck 



(Al) 



(A2) 



chain length L — An + 2, for which the ground state is 
non-degenerate with an odd number of particles, N — 
2n + 1. The Hamiltonian then reads 



3 3 



(A5) 



where e k = J cos k. The Bogoliubov transformation (k > 
0) 

c fe = cos -d k 7 fc + sin d k j k _ n , (A6) 
Cfe-7r = - sin $ k 7 fc + cos i9 fc 7 fc _ w , (A7) 



diagonalizes the Hamiltonian ( A5 ) , provided that 
tan2i?/j = h/e k . Thus we get 



where the single-particle spectrum is given by 



E k = sign(e fe )\/4 + h2 ■ 
The local moment rrij — (Sj) is found to be 

(-1)^ \- 1 
3 L ^ \E k \ 

k>0 1 1 



(A8) 



(A9) 



(A10) 



For L — >■ oo we replace the k sum by an integral and 
obtain Eq. Q. 



For even N, the undesirable sign of the boundary terms 
can be distributed in a homogeneous way over the entire 
ring using the transformation Cj — > 
tt/L. The transformed Hamiltonian 



e lvJ Cj. where <p — 



H=J 2 E( e< " c 5 c i+i + h.c.) - h J2 cfo (A3) 

3=1 3=1 



can then again be diagonalized by Eq. (Al), leading to 



H = [J cos (k + (p) - h] c^Cfc 
k 

= ^ (J cos k — h) c k c k , 



(A4) 



where in Y^' k the k values arc shifted with respect to those 
in J2 k ; they assume values ^r(v + |), —L/2 < v < L/2. 



2. Staggered field (L even) 

For a staggered transverse field, hj = (—l) J h, j = 
1, L, and an even number of sites the Hamiltonian has 
full particle-hole symmetry. For simplicity we choose a 



Appendix B: Perturbation theory for the step model 



We consider the step model, hj 



-h to the left and 



hj = h to the right. For simplicity we choose L = 2M 
with M odd. The ground state then has N = M parti- 
cles and the hopping term in Eq. ^ has the conventional 
tight-binding form. In two limiting cases, \h\ -C J and 
\h\ 3> J, we can calculate interface effects perturbatively. 
In these two limits, the Hamiltonian ^ is split in differ- 
ent ways into a dominant part Hq and a perturbation H' . 
For \h\ <C J we make the obvious choice 

j L 

H Q = - ^(ctcj+i + h.c.) = £ k c\c k , 



3=1 
M 



H' 



h ^2( c } c i ~ c l+i-j c L+i-j) 

3 = 1 



J2~ h 

k,k> 



k-k' c\c k ' , 



(Bl) 



where 



E k = J cos k , 

h„ = -ih(l - cos qM) . - 

sm(g/2) 



(B2) 
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while for \h\ » J we just interchange the two terms. 
Therefore the unperturbed ground states |$q) are, re- 
spectively, the half-filled band (levels with \k\ < § occu- 
pied, the others not) and the half-filled chain (left sites 
occupied, right sites empty). Using the fundamental re- 
lation of many-body perturbation theory, we can obtain 
the local moment rrij — (c\cj) 



^ from 



(cU> = (^ \T{c] Cj C/(oo,-oo)}|$ )c 



where 



U(oo, 



exp 



+00 



dt H'(t) 



(B3) 



(B4) 



T is the time-ordering operator and only connected dia- 
grams contribute. 

We discuss first the case of weak fields, \h\ <C J. In 
the unperturbed ground state the magnetic moments rrij 
vanish. The first-order correction is found to be 



rn 



(i) 



■5E 



L 2 



k,k' 



^^h^, 1 ^^, (B5) 
£k — £fc' 



where is the Fermi distribution of the unperturbed 
ground state. In the thermodynamic limit the sums 
are replaced by integrals. Introducing the variables 
Q = i(fc + k') and q — k — k' , we can perform one in- 
tegration and obtain for the region close to the interface 
located between sites M and M + 1 



The remaining integral is readily computed numerically. 
One obtains a step-like contribution corresponding to av- 
erage positive and negative magnetizations to the left 
and right of the interface, respectively, plus an alternat- 
ing part (wave vector ir) which decays like l/\£\ away 
from the interface. For sites close to the other interface 
(loca ted between the connected chain ends) the pattern 
(B6) is simply multiplied by a minus sign. This per- 



turbative formula agrees well with numerical results for 
h < 0.2 J. For larger values of h the dominant wave vec- 
tor is no longer 7r, but tt/ arccos (h/ J), as explained in 
Section IV. 

We now discuss the other limit, h 3> J, where the 
hopping term is the perturbation. In the unperturbed 
ground state the sites to the left of the interface are oc- 
cupied while those to the right are empty, corresponding 
to magnetic moments +5 and —5, respectively. It is 
easy to see that the first order correction vanishes. Up 
to second order we find 



m M +e 



1 

2' 

2 V 4/1 , 



e<-i, 
£ = o, 



-i + (jLy £ = i 
-I, e>2. 



(B7) 



mM+i 



(-iY 



dq 



cosq(£ — \) 



=2 9 



log(tan|). 



(B6) 
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